Strong Interactions of Single Atoms and Photons near a Dielectric Boundary 
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Modern research in optical physics has achieved quantum control of strong interactions between 
a single atom and one photon within the setting of cavity quantum electrodynamics (cQED) 1 . 
However, to move beyond current proof-of-principle experiments involving one or two conventional 
optical cavities to more complex scalable systems that employ N 1 microscopic resonator^ 
requires the localization of individual atoms on distance scales < lOOnm from a resonator's surface. 
In this regime an atom can be strongly coupled to a single intracavity photorP while at the same 
time experiencing significant radiative interactions with the dielectric boundaries of the resonator^. 
Here, we report an initial step into this new regime of cQED by way of real-time detection and 
high-bandwidth feedback to select and monitor single Cesium atoms localized ~ 100 nm from the 
surface of a micro-toroidal optical resonator. We employ strong radiative interactions of atom and 
cavity field to probe atomic motion through the evanescent field of the resonator. Direct temporal 
and spectral measurements reveal both the significant role of Casimir-Polder attraction^ and the 
manifestly quantum nature of the atom-cavity dynamics. Our work sets the stage for trapping atoms 
near micro- and nano-scopic optical resonators for applications in quantum information science, 
including the creation of scalable quantum networks composed of many atom-cavity systems that 
coherently interact via coherent exchanges of single photons^. 



The proximity of dielectric boundaries fundamentally 
alters atomic radiative processes as compared to quan- 
tum electrodynamics in free space. For example, free- 
space Lamb shifts and Einstein- A coefficients (i.e., level 
positions and decay rates) are modified for atom-surface 
distances comparable to the relevant transition wave- 
lengths, as considered in the pioneering analyses of 
Casimir and Polder 5 and of PurcelP in the late 1940s. 
Seminal experiments in the 1970s investigated radiative 
decay for organic dye molecules near a metal mirror^ 
and were followed in the 1980s by landmark observations 
of the inhibition of spontaneous emission for a trapped 
electrorP and an atom in a waveguide 9 . The ensuing 
years have witnessed the development of cavity quantum 
electrodynamics (cQED) in this perturbative regime of 
boundary-modified linewidths and level shift s^^J, with 
applications ranging from measurements of fundamental 
constants^ to the development of novel semiconductor 
device^. 

With increased interaction strength, a non-perturbative 
regime of cQED becomes possible and is character- 
ized not by irreversible decay but rather by the cyclic, 
reversible exchange of excitation between atom and 
photon 14 . The experimental quest for strong atom- 
photon coupling had its initial success in 1985 in the mi- 
crowave regime with the realization of the micromasei^, 
with strong coupling in the optical domain achieved some 
years latei™. By now the coherent control of atomic ra- 
diative dynamics has become possible on a photon-by- 
photon basi J 1|17 i Strong coupling has also been demon- 
strated for a wide class of physical systems^ beyond 
single atoms, including quantum dots coupled to mi- 
cropillars and photonic bandgap cavities^ and Cooper- 
pairs interacting with superconducting resonator J 2Q | 21 l 



This non-perturbative regime of cQED with strong light- 
matter interactions mediated by single photons has led to 
new scientific capabilities, ranging from a laser that oper- 
ates with one-and-the-same atoirP^ to the deterministic 
generation of entangled photon pairs^ to a two-qubit su- 
perconducting quantum processoiJ^. 

To a large degree, advances in the perturbative and 
non-perturbative regimes of cQED have been made inde- 
pendently. For example, for one atom localized near the 
center of a Fabry-Perot cavity with volume (7) 3 ~ (10 
/im) 3 , the coherent coupling g to an optical resonance 
can be large compared to radiative decay character- 
ized by the Einstein- A coefficient and cavity loss rate k, 
namely g ^> (70, k) where 70 = -A/2, placing the system 
in the regime of strong, non-perturbative atom-photon 
coupling^. Nevertheless, corrections to the atomic Lamb 
shift and Einstein- A. coefficient arising from surface in- 
teractions with the cavity boundaries remain small (e.g., 
5 A I 'A ~ 10 -5 ). However, many applications in Quan- 
tum Information Science could benefit from strong atom- 
photon interactions with micro- and nano-scopic optical 
resonator d 25 " 28 [ Atomic localization on a sub-wavelength 
scale near a resonator's surface is then required, with as- 
pects of both perturbative and non-perturbative cQED 
necessarily coming into play. 

In this manuscript we investigate such a regime for 
single Cesium atoms radiatively coupled to a high-Q mi- 
crotoroidal cavity 3 -^ and localized near the resonator's 
dielectric surface. As illustrated in Fig. [T^i, cold Cesium 
atoms are released from an optical dipole-force trap and 
randomly fall past the microtoroid. A real-time detection 
scheme based upon strong radiative interactions between 
one atom and the evanescent field of the cavity selects 
atomic trajectories localized within d < 300 nm from 
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FIG. 1. Radiative interactions and optical potentials 
for an atom near the surface of a toroidal resonator. 

a, Simple overview of the experiment showing a cloud of cold 
Cesium atoms that is released with then a few atoms falling 
near a microtoroidal resonator. Light in a tapered optical 
fiber excites the resonator with input power P[ n at frequency 
cj p , leading to transmitted and reflected outputs Pt,Pr- b, 
Cross section of the microtoroid at = showing the coher- 
ent coupling coefficient \g (f) — g(p,z,<fi)\ for a TE polarized 
whispering-gallery mode. The microtoroid has principal and 
minor diameters (D p ,D m ) = (24,3) /im, respectively, c, (i) 
Coherent coupling \g(d, z, (j>)\ for the external evanescent field 
as a function of distance d — p — D p /2 from the toroid's sur- 
face for (z,(f>) = (0,0). (ii) The effective dipole potentials Ud 
for resonant uj p = uoi -' , red uj p < cji ^ and blue u p > free- 
space detunings of the probe Pi n . The Casimir-Polder surface 
potential U B for the ground state of atomic Cs is also shown, 
(iii) The atomic decay rate 7(d) as a function of distance d 
from the toroid's surface for TE (711) and TM (7^) modes. 
All rates in this figure are scaled to the decay rate in free 
space for the amplitude of the Cs 6P3/2 — > 66*1/2 transition, 
7o/27r = 2.6 MHz. The approximate distance scale probed in 
our experiment is < d < 300 nm. 



the resonator's surface, with a large fraction of atoms 
passing below 100 nm and crashing into the surface. 
On this scale, the atom's coherent interaction with the 
cavity field is characterized by strong, non-perturbative 
coupling [Fig. [TJd, []J(i)], which we demonstrate by di- 
rect measurements of so-called "vacuum-Rabi" spectra 
for light transmitted and reflected by the atom-cavity 
system, as well as by observations of photon antibunch- 
ing for the transmitted light. On the other hand, the 
atom's motion and level structure are significantly in- 
fluenced by the (perturbative) Casimir-Polder potential 
from the surface's proximity [Fig. [Ij^ii)], which we infer 
from measurements of the time dependence of the cavity 



transmission during an atomic transit event, as well as 
from modifications of the spectra recorded for the trans- 
mitted and reflected fields. These observations are in 
reasonable agreement with a theoretical model that we 
have implemented by Monte- Carlo simulation and which 
gives insight into the underlying atomic dynamics, as de- 
tailed in Ref. [29l 

For the identification of atoms near the surface of the 
microtoroid in the regime shown in Fig. [TJ^, we rely on the 
strong interaction of atom and cavity field to modify the 
light transmitted by the cavity. Specifically, because the 
atom-cavity coupling coefficient g (r(t)) depends upon 
the atomic trajectory r(t), we can select single atoms lo- 
calized in the cavity mode by demanding a minimum 
criterion for the change in cavity transmission due to 
the atomic trajectory. Our scheme for single- atom d etec- 
tion is similar to that used in previous workP^^, but 
with significant modifications. Namely, by implementing 
fast digital logic, we achieve reliable real-time identifica- 
tion of atomic transit events in times as short as 250 ns 
from the photoelectric counts from the transmitted power 
.Pt(£)- Given the identification of a localized atom, the 
control logic switches the power P{ n and frequency u p of 
the probe input within ~ 100 ns and records subsequent 
photoelectric counts for the transmitted Pt(^) an d re- 
flected Pr(£) outputs from the cavity. These records of 
photoelectric counts form the basis for our analysis that 
follows, with further details presented in the Appendix 
and Supplementary Information (SI). 

To address experimentally the question of the distance 
scale for the recorded atom transit events, we first ex- 
amine the time dependence of the cavity transmission 
T(t) immediately following a trigger heralding the ar- 
rival of an atom into the cavity mode. Figure [2^i shows 
T(t) = Pr(£)/Pin for the case of resonant excitation, 
namely A pa = u p - cj a 0) and A ca = u c - u; a 0) 0, 
where uj^ is the free-space atomic frequency for the 
6Si/2,F — 4 — >• 6P3/27P = 5 transition in atomic Cs 
and uo c is the resonant frequency of the toroidal cav- 
ity. Two characteristic decay times are evident, with the 
background subtracted transmission Ts(t) = T(t) — B 
fitted well by the sum of an exponential (ex e -t / 5£/ ) 
with Sti = 0.78 ± 0.02 /is and a Gaussian (oc e -(V^n) 2 ) 
with Stu = 3.75 ± 0.09 /is. Here, the background level 
B = T(t ^> Stiji) is determined from the cavity trans- 
mission for times long compared to the duration of the 
transit event. 

The time constants 8t\, Stu can be associated with dis- 
tance scales di, du by way of the average velocity v with 
which atoms arrive at the toroid's mode following release 
from the optical trap, namely v ~ 0.17 m/s, leading to 
d\ ~ 130 nm and d\\ ~ 640 nm. For comparison, the scale 
length for g(d) in the radial direction is A = l/k$ = 136 
nm (Fig. [^(i)), while in the vertical direction, the vari- 
ation of g(z) is approximately Gaussian (oc e~^ z l W{ ^ ) 
with waist wq — 590 nm (Fig. |T]o) . The comparisons 
d\ ~ X and d\\ ~ wq suggest that the short-lived compo- 
nent Sti in Fig. [2^i arises from atomic trajectories that 
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FIG. 2. Observation (a) and simulation (b-e) of atomic 
transits within the evanescent field of the micro- 
toroidal resonator for A ca = A pa = 0. a, Observed 
cavity transmission Tb (£) versus time t following a triggering 
event at t = 0, with approximately 5 x 10 4 triggered tran- 
sits included. The data are fit to the sum of an exponen- 
tial (I) and a Gaussian (II) (red curve), with time constants 
8tj = 0.78 ± 0.02 fjs and 5t u = 3.75 ± 0.09 fis, with each 
component shown by the dotted lines, b, Simulation result 
for the cavity transmission (£) versus time t (points) from 
an ensemble of triggered trajectories. The red curve is a fit 
to the sum of an exponential and Gaussian with time con- 
stants St[ s ^ = 0.69 /is, St[f = 4.0 /is while the dotted lines 
represent the individual fit components, c-e Probability den- 
sities pi (d), pi (g), pi (5 a) for the distance d, coupling g, and 
transition frequency shift 5& = uo a (d) — from the same 
simulation set as for (b). {d, g, 5 a } are averaged over the first 
500 ns following the trigger. For these results, the trajectories 
are divided into two classes, i — {I, II} corresponding to the 
two time constants 8t[ s ^ (blue shaded curve) and 6t[f (red 
shaded curve) in (b). 



are deflected from their otherwise vertical fall to largely 
radial paths of length A that terminate at the dielec- 
tric surface. Similarly, the longer-lived component Stn 
is associated with trajectories that pass along z without 
significant radial motion toward the surface of the toroid. 

Of course an atom near the surface will not move 
with constant velocity but will be accelerated by inter- 
actions with surface potentials and the cavity field itself. 
To reach a quantitative understanding of the external, 
center-of-mass motion and the internal, atomic dipole- 
cavity field coupling, we have implemented a numeri- 
cal simulation that incorporates both perturbative and 
non-perturbative aspects of the radiative interaction of 
the atom and micro-toroid. Our Monte-Carlo simulation 
draws random initial trajectories for atoms falling from 
a thermal cloud and implements a stochastic process for 
photoelectric detection to emulate our real-time detec- 
tion technique. The model includes Casimir-Polder and 
dipole forces from the potentials U s (r),Ud(r) shown in 
Fig. Ilb(ii), atomic level shifts (and hence detunings) from 
?7 s (fj7and boundary-modified decay 7||(d) shown in Fig. 
[jj(iii). The non-perturbative interaction of atom and 
cavity field is based upon the analytic results in RefP^. 
Details of the simulation appear in Ref. [29] and the SI. 

Results from this analysis are presented in Fig. ^[b-e) . 
In agreement with the observations in Fig. |2^i, Tb (t) from 
the simulation in Fig. [2}d exhibits two time scales and is 
fit well by the sum of an exponential with time constant 
St[ s ^ = 0.69 [is and a Gaussian with 1/e width Stjj = 4.0 
/is. Atomic trajectories associated with the shorter time 

(s) 

Stj have distances peaked around d\ ~ 100 nm and ter- 
minate with crashes into the surface of the toroid [pi(d) 
in Fig. [2^], but exhibit large coupling gi/2ir ~ 40 MHz 
\pi(d) i n Fig. and large surface- induced shifts of 

the atomic transition frequency S aj i/27T > 10 MHz for 
d < 60 nm [pi(S a ) in Fig. [2^]. By contrast, the com- 
plementary set of trajectories for the longer time St^j 
pass roughly vertically through the cavity mode along 
z (hence their Gaussian time dependence). Relative to 
the St^ set, these trajectories exhibit larger distances 
d\\ ~ 250 nm, smaller coupling gijl-K ~ 20 MHz, and 
surface-induced shifts of the atomic transition frequency 
<W < 7o \pn(d) : pii(g) : pii(S a ) in Figs. [2|c-e), respec- 
tively] . 

In Fig. [3] we investigate temporal dynamics for the 
cavity transmission T(t) but now with non-zero detuning 
between the atom and probe field, A ca /27r = A pa /27r = 
±40 MHz (Fig. 3a). Since uj p ^00^ dipole forces from 
coherent excitation of the intracavity field should induce 
an asymmetry for T(t) for red and blue detuning, with 
faster decay for red detuning (uj p < uj a ) due to the com- 
bined effect of the attractive potentials U s (f) and U&(r) 
shown in Fig. []J(ii). The data are fit well by decaying ex- 
ponentials (ex e~ l l bti ) added to Gaussian components in 
analogy with the earlier analysis, with exponential time 
constants 5t re( ± — 0.11 ± 0.01 and ^tbiue = 0-53 ± 0.03 
/is. Also plotted are simulations of the experiment for 
freely falling atoms with U s = Ud = 0, atoms under the 
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FIG. 3. Dynamics and trajectories for strongly coupled atoms moving in surface and dipole potentials {U s , Ud}- 
a, Transmission T(t) for A ca /27r = —40 MHz (left) and +40 MHz (right) measured after an atom trigger at t = 0. In each 
panel, the circles are data for 2 x 10 3 trigger events; the lines are simulations of T(t) for the full model (blue), for U s = 
(magenta), and for U s = Ud = (g reen) . b ? Representative atomic trajectories projected onto the p — z plane for simulations 
in panel (a), with the TE mode intensity plotted on a gray scale. The upper panels are for A C8l /2tv = —40 MHz while the 
lower panels are for A ca /27r = +40 MHz. The color bars at the top of the panels match the colors of the curves in (a). For 
each panel, orange lines are untriggered trajectories, while triggered trajectories are represented by blue lines which turn red 
after a trigger at t = 0. c, Simulations showing trajectories from a full 3D simulation with U B , Ud, as well as a two-color dipole 
potential (FORT) triggered "on" by atom detection at t = 0. A ca /27r = +40 MHz in correspondence to (a), (b). Blue lines 
represent falling atoms with the FORT beams "off' (t < 0), while red lines are trajectories after the FORT is triggered "on" 
and an atom begins to orbit the toroid. To illustrate the timescale, the trajectories are colored pink for t > 50 ps. 



influence of only dipole forces, U s = 0, and a full simula- 
tion including both the dipole force and CP forces. For 
both red and blue detunings, the timescales from the 
simulation absent U s ,Ud are substantially longer than 
observed in experiment, while the full simulation with 
U s , Ud leads to time constants St^ d = 0.19 ± 0.02 fis and 

^mL = 0-59 ± 0.06 {is. These simulated time constants 
are qualitatively similar to observations, with quantita- 
tive differences attributed to simplifications inherent in 
the simulation model (see SI). In contrast with the case 
<^p = in Fig- |2j the Gaussian component of these tem- 
poral decays is minimal because of the difference in scale 
lengths of the CP potential U s and effective dipole poten- 
tial Ud(A ca = ±157o), which become comparable to 70 
for distances d < {65, 200} nm, respectively (Fig. [^(ii)). 
Long range dipole forces, which are largely absent for 
Fig. [2j dominate the trajectory dynamics of Fig. [3] and 
consequently vertically falling long-lived Gaussian trajec- 



tories do not significantly contribute. 

To illustrate the underlying atomic motion, Fig. ^jp dis- 
plays atomic trajectories projected onto the p — z plane 
for each simulation in Fig. [3^i. Each panel displays a 
representative sample of untriggered and triggered tra- 
jectories. For red detuning, introducing dipole forces and 
CP forces leads to every triggered atom crashing into the 
toroid surface, explaining the short decay 5t re d. The blue 
detuned case is more complicated, with both attractive 
CP forces and the repulsive dipole force reducing the time 
the atom is in the mode; CP forces pull nearby atoms into 
the surface while the dipole force repels other atoms out 
of the mode. 

As shown in Fig. [3J3 and discussed in the Appendix, 
we have augmented our numerical simulation to include 
a dipole force optical trap J7 t (FORT) formed by the 
toroid's evanescent fielcpMM] ^ n addition to the poten- 
tials U s , Ud- The trapping potential U t is triggered on 
by the same criteria as for Figs. [3^i, b, with then a signif- 
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icant fraction of triggered atoms bound in orbit around 
the toroid for durations surpassing 50 /is. 

The measurements in Figs. [2j [3] rely upon strong in- 
teractions of one atom and photon for initial atomic lo- 
calization within the cavity mode and for measurements 
of the subsequent motion by way of T(t). To establish 
directly the non-perturbative coupling of atom and cav- 
ity field, we next turn to measurements of transmission 



(a) 



(b) 



T(cjp) = PT(^p)/^in and reflection R(u 



PrH)/Pit 



spectra as functions of the frequency uo p of the inci- 
dent probe field P in . Probe spectra {T(cj p ), R(u p )} are 
recorded following the detection of a single-atom event 
with co> p = uo c for a fixed detuning A ca between atom 
and cavity to optimize sensitivity for an intracavity atom 
(i.e., A ca = A pa ). With an atom thereby localized in 
the cavity mode, fast control logic and feedback switch 
the probe power P in to some fiducial level for a given 
spectrum and the probe frequency uj p to a relevant de- 
tuning A pa 7^ A ca for measurements of {T(cj p ), R(uj p )}. 
Because falling atoms remain in the evanescent field for 
only a few /is, the spectra are built up over thousands of 
transit detections and consequently represent an ensem- 
ble average over triggered atom trajectories. 

Strong radiative coupling of an atom and a mi- 
crotoroidal resonator is described by an extension of 
the Jaynes-Cummings 14 Hamiltonian 3 (see SI). Our 
whispering-gallery resonator supports two counter- 
propagating traveling- wave modes {a, b} that are cou- 
pled by internal scattering at a rate /i, leading to cav- 
ity eigenmodes {A, B} that are linear superpositions of 
{a, b}. The interaction of the {A, B} modes with an 
intracavity atom is characterized by coherent coupling 
g(r), with the resulting atom-cavity eigenvalues {A^} 
shown in Fig. [4^i for the single-excitation eigenstates. 
For large detuning |A ca | ^> g, there is one atom-like 
and two cavity-like (for the {A, B} modes) eigenvalues. 
For 

^ g~> there is the familiar anti-crossing between 
the imaginary parts of two dressed-state eigenvalues X± 
with splitting AA ± = Im (A+ - A_) « + 4 # 2 for 

g ^> {ft, ft, 7||}, while the third cavity-like eigenvalue Ao 
remains uncoupled to the atom. This dressed-state eigen- 
structure, along with the dissipative rates J\\(d) and n 
for atom and cavity, determine the system's spectral re- 
sponse {T(cjp), R(u p )}. 

Using a simple model with atoms falling vertically 
through the evanescent field of Fig. lb with {Z7 S , Ud} — 
(SI), we construct a probability distribution Pf a ii(#) of 
coupling constants for atom detection, with probe spec- 
tra {T(cjp), R(u p )} then obtained by averaging spectra 
for fixed g over the distribution pf a n(g) (Fig. 4b). Al- 
though the full eigenstructure from Fig. 4a cannot be 
resolved due to the 'smearing' from PM\(g) even with 

#max > {«,7||}j tne splitting Acj^ aks between A_ and 
Ao (shifted by its proximity to the unresolved A+) is 
resolved and approximates, though underestimates, the 
eigenvalue splitting AA± (i.e., Acjp^ ks /27r = 110 MHz, 
while AA ± /2tt = 130 MHz). 

Figure]^, d show measured spectra for both the bare- 
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FIG. 4. Transmission T(uj p ) and reflection R(u p ) spec- 
tra for single atoms coupled to a microtoroidal res- 
onator, a, cQED eigenvalues A±,o for {h,g}/27v = {10,40} 
MHz as a function of atom-cavity detuning A ca . The dashed 
lines indicate the detunings for the spectra in the follow- 
ing panels, b, T(cj p ) for A csl /2tv = +60 MHz from a sim- 
ple average for falling atoms over the distribution pfaii(g) 
(inset) absent cavity and surfaces forces. Acj pe aks is found 
from the peaks indicated by arrows, c, d, Experimental re- 
flection Ri(A pa ) and transmission Ti(A pa ) spectra for the 
empty cavity i = NA (red) and with atoms i = A (blue) 
with the peaks used for Acj e x P indicated. Curves are re- 
sults of the full Monte Carlo simulation, e Difference spectra 
AR = R A (A pa ) -i?N A (A pa ) and AT = T A (A pa ) - T N A(A pa ) 
for A ca /27r = +60 (i,ii), +40 (iii,iv), -40 MHz (v,vi). Green 
lines are simulation results for U s = Ud = 0, while blue lines 
are from the complete simulation. 



cavity with no atoms (NA), -Rna (A pa ) and Tna (Apa), 
and with triggered single atoms (A), Ra (A pa ) and 
Ta (A pa ), for A ca /27r = 60 MHz. The splitting 
Au; ex p/27r ~ 95 ± 5 MHz between the prominent cavity 
peak and the dressed state feature can be read directly 
from both Ta and Ra- Taking Aco> p ^ s as a lower es- 
timate for the average eigenvalue splitting AA± yields 
an average coupling g/2ii « 37 ± 3 MHz. This average 
coupling indicates that strong coupling is achieved on 
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average, with g > (ft, 70), where (ft,7o)/27r = (21,2.6) 
MHz. 

Quantitative differences between the simple model in 
Fig. [4)3 and the experimental spectra in Fig. [4J3, d yield 
information about effects beyond the simple model, in- 
cluding perturbative surface interactions not included in 
the standard Jaynes-Cummings treatmentP^. In par- 
ticular, the feature at A pa /27r = —30 MHz in both Ra 
and Ta is significantly reduced in magnitude from the 
spectrum predicted by Pf a ii(#), which as discussed below, 
results from the effects of {[/ s , U^} on the atomic trajec- 
tories and internal levels as in Figs. [2j [3j 

Measurements of the difference spectra taken with and 
without atom transit events, Ra — Rna and Ta — Tna 
are shown in Figs. [4^ for cavity detunings A ca /27r = 
60,40, —40 MHz. Again, the simple prescription of read- 
ing Acjp ea ks directly from the splitting of the low and high 
frequency peaks together with the expression for A\± 
leads to an estimate of the average coupling g/2it & 35 ±5 
MHz that is consistent across the six spectra displayed. 

For comparison to the measured spectra, the full curves 
in Figs. [4^ are from our Monte Carlo simulation for 
AR = iJ(Ap a ) - iV>o(A pa ) and AT = T(A pa ) - 
T^o(Ap a ). Calculated spectra are shown both for the 
full model and with all forces removed. Agreement with 
the full model is achieved for the choice g m ax/27r ~ 100 
MHz, which is somewhat less than the value of g m ax/27r = 
140 MHz expected for the fundamental TE mode near 
A = 852 nm estimated from a finite element calculation 
(Fig. [TJd). The difference is attributable to imprecise 
knowledge of the toroid geometry and mode. Except for 
the relevant detunings and measured cavity decay rates, 
the same parameters are used for each spectrum simula- 
tion; specifically, (^ max , 70) /2tt = (100,2.6) MHz. Note 
that apart from the adjustment of ^ max , all parameters 
in the simulation are estimated from measurements of 
our system or, in the case of dipole forces, surface forces, 
and surface level shifts, are taken from theoretical and 
experimental results in the literature (SI). 

For each of the spectra in Fig. [4^, removing the 
Casimir-Polder and dipole forces (i.e., setting U s ,Ud to 
zero) leads to increased deviations from the measured 
spectra relative to the full simulation, which describes the 
measurements reasonably well. The most significant ef- 
fect of U a is seen for a red-detuned cavity (A ca /27r = —40 
MHz) where significant spectral features not readily ob- 
served in the data appear for U 8 = 0. In combination 
with the temporal analysis in Figs. [2j [3J the cQED 
spectra illustrate the necessity of including perturbative 
surface interactions for understanding atomic dynamics 
near the resonator. The model uses a distance depen- 
dent atomic decay rate 7||(d) for the linearly polarized 
TE mode, but the differences between 7y and 70 are too 
small to be observed in the data. Despite the overall 
consistency achieved with the full simulation, systematic 
disagreements between data and model suggest that fur- 
ther analytical progress is required, including better (in- 
dependent) knowledge of the toroid geometry, as well as 
solving the full master equation to account for the multi- 



level structure of the Cs atorrP^. 

To confirm the quantum nature of the atom-cavity in- 
teraction near the surface of the toroid, we present in 
Fig. [5] measurements of photon statistics for the trans- 
mitted field Pt for A ca = 0. Photon statistics are in- 
ferred from the time records Ci^(U) of photoelectric 
counts in time bins U for two detectors D\ , D2 follow- 
ing an atomic trigger event at t = 0. Fig. [5] shows the 
average cross-correlation Ci2(r) = Y^i(Ci(U)C2(U + t)) 
for < ti < 8 /is as well as the cross-correlation of the av- 
erage counts Ci 2 (t) = ^i{Ci(ti)){C2(ti + T)), where the 
angled brackets represent sums over the ensemble of trig- 
gers. The photon process is super-Poissonian, indicated 
by C\2{r) > Ci2(r) for all r, due to large fluctuations in 
transmitted intensity from one atom to the next, which 
presumably arise from variations in atomic position r 
and hence coupling g(r) near the toroid's surface (in- 
set (i)). Even in the face of these large fluctuations, the 
non-classical character of the atom-cavity interaction sur- 
vives, as is evident from the short-time dynamics shown 
in Fig. [5j where 612(0) < C\2(r) exhibits photon anti- 
bunching. From the minimum at r = 0, C\2(r) regress 
to its peak with characteristic half-width of 6 ns. 

In contrast to the case C12 —> realized in microtoroids 
in the bad-cavity limit for a photon turnstile 30 and a pho- 
ton router 3 ^, here 612(0) is 0.55 of its maximum value. 
This behavior is a result of coherent dressed-state dy- 
namics for A ca = rather than a projective measure- 
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FIG. 5. Photon statistics for atoms with A ca = 0, A pa = 

0. Cross-correlation Ci2(t) (blue circles) computed from the 
records of photoelectric counts at detectors Di,D 2 from the 
forward flux Pt from a sum over many atom trajectories show- 
ing photon antibunching around r = 0, with C±2(t) obtained 
from the product of averages of the recorded counts at each 
detector for comparison (black circles). The red curve is a 
calculation for the two-time second-order correlation function 
averaged over the distribution p(g) from the full simulation 
and scaled by a single parameter to match Ci2(r) at r = ±40 
ns. (i) Expanded view of Ci2(r) and Ci2(r) over full range 
of r, with the long decay related to the typical atom transit 
time. 
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ment of the atom as in Refs. [30] and [3TJ Averaging 
a time-dependent calculation for the coincidence count 
rate for a fixed value of g over the distribution of g ob- 
tained from our trajectory simulations results in reason- 
able agreement with our measurements (red curve in Fig. 
[5]). Except for the assumptions inherent in our simula- 
tion, the only free parameter in this calculation is the 
amplitude which is scaled to match the data. Our model 
predicts both 612(0) =^ and the regression timescale 
near r = 0, which supports its effectiveness in describing 
the quantum behavior of the atom near the surface of the 
toroidal resonator. 

By exploiting real-time triggering of single atoms, our 
experiment has realized a system where an atom's dy- 
namics are governed by both its strong, single-photon 
interactions with the resonator's field and perturbative 
forces on classical atomic motion and internal level struc- 
ture from proximity to the resonator's surface. Enter- 
ing this regime opens the door for quantitative study of 
dynamical Casimir-Polder forces in the strong-coupling 



limit 36 37 , which will require trapping atoms at short dis- 
tance scales for extended interrogation times. In contrast 
to the standing- wave structure of a Fabry- Perot cavit}P^, 
mirotoroidal resonators offer the tantalizing possibility 
of radiall y conf ining an atom in a circular orbit around 
the toroicP^MJ, w ith initial results from our simulation 
shown in Fig. [3J3. In correspondence to the development 
of cQED to reach the regime of strong coupling with one 
trapped atom in a Fabry- Perot resonator 38 , the advances 
described here offer an important advance toward trap- 
ping and cooling of a single atom near the surfaces of 
micro- and nano-scopic optical resonators. 
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Appendix 

Silica microtoroidal resonators with principal diame- 
ter D p ~ 24 /im and minor diameter D m ~ 3 /im 
are fabricated using standard methods 26 and mounted 
in an ultra-high vacuum chamber with pressure below 
10 -9 torr. The cavity resonance frequency uj c is tuned 
using the silicon substrate temperature to be near the 
6Si/2,F = 4 — » 6P 3 / 2 ,P' = 5 transition of Cs at fre- 
quency Cesium atoms are cooled below 10 fiK in 
a magneto-optical trap and transported by an optical 
conveyor belt ~ 800 /im above the toroid!^ The mi- 
crotoroid is excited by a tapered fiber carrying an in- 
put power Pi n ~ 4 pW at frequency u; p , which is cou- 
pled to a whispering-gallery mode of the toroid using 
piezoelectric positioners. The cavity internal losses ^i, 
fiber-taper coupling ft ex , and intermodal coupling h are 
{^i, ^ex, h}/27r w {8,10,13} MHz, while the maximum 
atom-cavity coupling rate at the surface of the toroid is 
inferred to be g m ax/27r ~ 100 MHz (i.e., a maximum 

single-photon Rabi frequency £lmL/27r = 2g max /27r « 
200 MHz). The transmitted (P T ) and reflected (P R ) 
powers from the coupling fiber are each directed to two 



single photon counting modules (SPCM), {Di,D 2 } and 
{Di rj P>2r}, respectively. At critical coupling^, the input 
field in the fiber taper and the field coupling out of the 
cavity destructively interfere in the forward transmission 
direction, so that Pr < 0.01Pi n for co p = u c . An atom 
coherently coupled to the cavity at rate g(r) disturbs this 
critical coupling condition and the photodetector counts 
at P>i,P>2 increase. Single photoelectric events within 
a running time window At tn are counted and compared 
with a threshold number Cth by a field programmable 
gate array (FPGA) operating at 40 MHz. For the data 
shown in Figs. 3 and 4, the parameters At t h = 750 ns 
and Cth = 5 are chosen to give a false detection rate less 
than 1% and an average trigger time as early as possible 
during the atom transits (SI). When the threshold crite- 
ria is met (defined as t = for each event), the FPGA 
sends a trigger pulse to a photon counting card to time 
stamp and record subsequent photodetections with 2-ns 
time resolution. 

At t = 0, the FPGA also triggers high-frequency op- 
tical modulators that switch the power and frequency of 
the probe input to the tapered fiber. For typical experi- 
ments, the probe flux is reduced to Pi n ~ 2 pW and the 
probe detuning A pa = uj p — uj^ set to a value within the 
range shown in Fig. [4j Including both optical and elec- 
trical delays, the optical switching is complete by t = 100 
ns. The photocount record on detector z, in a time 

interval < t < 8 /is following the trigger is recorded 
for a succession of N ^> 1 trigger events. For spectral 
measurements, the transmitted and reflected photocount 
records are averaged over a selected time window (typi- 
cally, 200 < t < 700 ns) and normalized by the photo- 
counts taken with large detuning, Pr (A pa ^> ft), to ob- 
tain the experimental transmission and reflection spectra, 
T (A pa ) and R(A pa ). Error bars for data are estimated 
assuming Poissonian counting statistics and are written 
as plus or minus one s.d. Fit results are quoted with 68% 
confidence intervals. 

Having validated our trajectory simulation with the 
measurements in Figs. [2][5j we have studied loading of 
falling atoms into a two-color ev anesce nt field far off- 
resonance trap (FORT) (Fig. ^jPtSEni A trapping po- 
tential Ut can be formed using a blue-detuned fundamen- 
tal mode and a higher order red-detuned mode^EHl p or 
our simulation, we use a red (blue)-detuned mode near 
898 nm (848 nm) with powers ~ 50 fiW each to give a 
trap depth of ~ 1.5 mK which is switched on at t = 
conditioned on a falling atom FPGA trigger. Despite the 
large kinetic energy of falling atoms and poor localiza- 
tion of the atoms relative to the trap minimum, approx- 
imately 20% of triggered atom trajectories are captured 
in the trap. Simulated trapping times exceed 50 /is, lim- 
ited not by heating from trapping light but by the ra- 
diation pressure from unbalanced traveling whispering- 
gallery modes. This radiation pressure leads to atom 
gallery orbits around the toroid (Fig. (3^). Exciting a red- 
detuned standing wave would provide three-dimensional 
trap confinement and increase the trap lifetime. 



